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Abstract. Stability properties are presented of dipole and quadrupole nonradial oscillation modes of model stars that experi- 
enced a late helium shell flash on their way to the white-dwarf cooling domain. The computed instability domains are compared 
with the observed hot variable central stars of planetary nebulae and the GW Vir pulsators. 
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1. Introduction 

The family of pre-white dwarfs (PWDs) referred to as PGl 159 
stars is spectroscopically defined by the dominance of He II, 
C IV, and O VI (and sometimes N V) lines and a strong H defi- 
ciency. Surface mass abundances of about 33% He, 50% C and 
17% O are reported as typical. Currently, 32 stars are members 
of the PGl 159 family, the effective temperatures range from 
about 170 000 to 80 000 K. Roughly half of the PGl 159 stars 
are embedded in a planetary nebula. 

Until recently, a hydrogen-deficient composition of the 
PGl 159 kind could not be reconciled with stellar evolution sce- 
narios. Only the introduction of appropriate overshooting and 
mixing procedures in modeling post-AGB stars undergoing a 
very late thermal pulse (VLTP) led to PWDs with H-deficient 
surface layers that are in agreement with PGl 159 spectroscopy 
llHerwigl200lL and references therein). 

Ten of the PGl 159 stars are pulsating variables (also bap- 
tized as GW Vir variables, after the variable-star designation 
of the prototype PGl 159-035; we will use both names syn- 
onymously throughout the text). Four of the ten variable stars 
have high surface gravities (logg > 7) and low luminosities, 
i.e. they are essentially on the corresponding white-dwarf cool- 
ing track and have hence evolved around the "evolutionary 
knee" (cf. Fig.^. The rest of the GW Vir stars has low logg 
(i.e. logg < 7), in four cases associated planetary nebulae are 
known. The observed periods range from about 5 to some 17 
minutes on the high-log g branch and they can exceed one hour 
in the low-log g PGl 159 variables. Based on these time-scales, 
the pulsations must be attributed to intermediate to high-order 
g modes. 
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The excitation physics of the PG1159-type oscillations is 
still a matter of debate; this, despite the relative simplicity of 
the microphysics in the destabilization region as compared with 
e.g. classical pulsating variables. Most of all, in the case of the 
PWD pulsations, the stability analyses are not discredited by 
the insufficiently understood role of convection which is usu- 
ally present when the /c-mechanism is at work. Nonetheless, 
no general agreement exists regarding in particular t he chemi 
cal compos i tion in the driving zone (e .g. Starrfiel^^t^ 
'Saio' '1996* 'Bradle v & Dziembowskil 11996.: .Gautschvl 
Cox 2003; Ouirion et alJ2004h . 

LSaiol ill 9961) and lGautsch\J l ll997h studied qualitatively the 
influence of the new generation of opacity data, specifically of 
the OPAL opacities, on the excitation of GW Vir pulsations. 
Due to the lack of evolutionary models with realistic histories 
that led them through the thermally pulsing AGB phase, the au- 
thors performed their explorative computations on model stars 
that evolved off the helium main sequence. The evolutionary 
phases before the knee differed significantly from post-AGB 
tracks. The post-knee epochs though, converged surprisingly 
well towards the loci of full evolutionary computations. The 
PG1159-type instability regions deduced from the simplified 
models agreed decently well with observations. The computed 
period separations, however, differed sufficient ly from the ob- 
served ones to arouse considerable criticism ( IO'Brienll2000l: 
■Cox,2003) . However, it seems to have largely gone unrecog- 
nized that usually the excitation physics is much more r obust 
than the asteroseismic signatures and that the studies of LSaiol 
(>1996..) and .Gautschv. L 1 997 .) served a useful purpose nonethe- 
less. 

In the following paper we return to the excitation physics 
of PGl 159 stars. We set out to show that the Saio and Gautschy 
results are robust, as it has also been shown recently in an in- 
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Fig. 1. Evolutionary track of the initially 2.7M0 star passing 
through the AGB, a very late thermal pulse (VLTP), and enter- 
ing the final cooling phase as a hydrogen-deficient white dwarf 
with O.59M0. Pulsation computations were performed at late 
epochs along the locus tracing the terminal departure from the 
red-giant domain. 



dependent approach by Ouirion et alJ ( l2004l) . This time, we go 
beyond the use of simplified evolutionary models or the use of 
envelope models to describe GW Vir pulsations. Numerical sta- 
bility analyses were performed on state-of-the-art model stars 
that evolved through a late helium shell flash after they evolved 
through the thermally pulsing AGB, starting originally from 
the main sequence. Hence, this paper closes the credibility gap 
left open by the earlier work. Section |2] introduces the stellar 
models and reiterates on the computational tools used for the 
stability analyses. Section|3presents the results of the stability 
computations for dipole and quadrupole modes; in particular 
the instability domains are described. The results are discussed 
and compared with earlier work in Sect. |4l the conclusions in 
Sect.|5]point out further steps that need to be taken. 

2. Models and Tools 

The model stars on which stability analyses were per- 
formed were extracted from the computations described in 
[AllJiaus et al . (2005); they foflowed an initially 2.7M0 star 
with solar composition from the zero-age main sequence, 
through the AGB and a VLTP, to the final cooling sequence 
as a hydrogen-deficient white dwarf. Except for the final 



white-dwarf cooling phase, the admixture of heavy elements 
in the star's envelope was assumed to be solar; the OPAL 
opacity tables that entered the computations are described in 
lAlthaus etai]J2005i) . The abundance evolution of 16 chemical 
elements was included via a time-dependent numerical scheme 
that simultaneously treated nuclear evolution and mixing pro- 
cesses due to convection and overshooting. Such a treatment is 
particularly important during the short-lived phase of the VLTP 
and the born-again episode for which the assumption of in- 
stantaneous mixing is inadequate. In the course of the VLTP 
phase, most of the residual hydrogen envelope is engulfed by 
the deep helium-flash induced convection zone and is com- 
pletely burned. The star is then forced to evolve rapidly back 
to the red-giant region and eventually into the domain of the 
central stars of planetary nebulae at high effective temperatures 
but this time as a hydrogen-deficient, quiescent helium-burning 
object. We note that the inclusion of overshooting below the 
helium-flash convection zone leads naturally to surface abun- 
dances in agreement with those observed in PGl 159 stars (see 
Fig.O. In particular, also the N abundance is compatible with 
the range of observed ones ( of the order of 0.01 in mass 
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pulsating PGl 159 stars (e.g. iDreizler & Heber..l998: .Werneil 
I2OOI). 

Four model-star sequences, with 0.53, 0.55, 0.59, and 
O.64M0, were analyzed in this paper. The O.59M0 se quence de- 
rived d irectly from the evolution computations of Alt haus et alJ 
('2005). Figure n displays the evolutionary track of the model 
with 2.7M0 on the zero-age main sequence and O.59M0 at 
the end on the white-dwarf cooling track. The stellar models 
with 0.53, 0.55, and O.64M0 were derived from the O.59M0 se- 
quence by changing the stellar mass to the appropriate value 
shortly after the termination of the born-again red-giant phase 
of the VLTP (i.e. at log Teff w 4.2). Upon restarting the evo- 
lution of the models with the abruptly changed mass, the ef- 
ficiency of the helium-burning shell did not match the model 
structure. For all three stellar masses, the transition phase dur- 
ing which nuclear burning and the internal structure re-adapted 
lasted to about log Tss = 5.1 on the locus of the last descent 
from the red-giant domain. 

Furthermore, the evolutionary phase along the white-dwarf 
cooling track (when logTeff < 5.0) included elemental diffu- 
sion which is important for the occurrence of a red edge of the 
pulsations as discussed in Sect. 13.21 

The nonradial stability analyses of dipole and quadrupole 
g modes were performed with the Riccati integration method 
of the linear no nadiabatic oscilla tion equations as used in the 
PGl 159 study o flGautschvl(ll997l) ; the approach was described 
for example in iGautschv et al.' ('19961 A gain, we retained the 
nuclear terms in the stability equations to detect potentially ex- 
istent e-driven instabilities. Convective energy transport is not 
important in the envelopes of the PGl 159 stars. Hence, the 
stability computations neglected convection effects completely. 
The Riccati method was described on various occasions in the 
past; nevertheless, elementary misjudgments continue to pre- 
vail (e.g. Cox 2003) and it seems necessary to emphasize once 
more that the Riccati approach is a shooting method. The spa- 
tial resolution that can be achieved in the pulsation compu- 
tations is not limited by the number of spatial gridpoints ap- 
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Fig. 2. Representative spatial abundance profiles for selected 
nuclear species for a O.59M0 model at log Teff = 5.15 and 
logL/Lo - 2.377. The outermost 10"'' M, are not shown as 
the compositions do not change there anymore; this applies 
to the evolutionary phases prior to logTeff = 5.0 along the 
white dwarf cooling track. Overshoot episodes during core he- 
lium burning and the AGB evolution left clear signatures in the 
abundance profiles: the step at -0.3 and the mountain-like fea- 
ture at around -1 of the abscissa. 

proximating the stellar background model. All that is required 
is that the stellar model must resolve all the relevant physical 
quantities sufficiently well. The physical quantities that must 
be accessible to the pulsation integration at arbitrary spatial co- 
ordinates are then interpolated in the stellar background via a 
monotonized cubic polynomial (Steffen 1992). The limitation 
of the radial mode order that was computed in this project (and 
any other before) was inflicted by the computing time that had 
to be invested per mode rather than by any numerical spatial 
resolution issue. 

In this project, the Brunt-Vaisala (BV) data had to be spa- 
tially smoothed. Very crisp composition steps in the stellar 
interior, be they imposed by elaborate treatments of convec- 
tive boundaries or artificially produced by the finite-difference 
methodology, cause spikes in the BV frequency. Shooting 
methods are particularly sensitive to spikes in the coefficients 
of the differential equations as they usually entail very short in- 
tegration steps to spatially resolve such rapid changes. Arguing 
that such discontinuities are somewhat smeared out in a real 
star, we smoothed the BV frequency numerically to iron out 
the spikes but maintaining at the same time its broader fea- 
tures. The simple algorithm that was invoked computed the BV 
frequency iteratively at any point in the star as the average over 
its nearest neighbors and consequently removed narrow spikes 
efficiently. 

3. Results 

This section describes the numerical results obtained from the 
nonadiabatic eigenanalyses of ^ = 1 and 2 modes, performed 
on the four model-star sequences with 0.53,0.55,0.59, and 
0.64Mo. 



Figure|3displays the instability boundaries on the log T^ff - 
period plane for t - \ (left panel) and £ - 2 modes (right 
panel). Eigenmodes with periods up to 1800 s were consid- 
ered. For the dipole modes, the observed domains which are 
occupied by GW Vir variables are overlaid for comparison. The 
ob served dipole-mo de period ranges were adopted from Tab. 1 
of iNagel & Werne'3 (i2004.) : the horizontal sizes of the boxes 
reflect the canonical 10 % uncertainty in the spectroscopic T^ff 
determinations. 

The log Teff - period diagram for dipole modes shows that 
the instability regions are separated into a low-log g (high lu- 
minosity) and a high-log § (low luminosity) domain for each 
mass. For all stellar masses considered, the tip of the evolu- 
tionary knee is pulsationally stable. The low-logg instabiUty 
domains are roughly compatible with observed periods in vari- 
able planetary-nebulae nuclei (PNNV). One exception is the 
star Lo4 with its short periods. The Lo4 observations could 
possibly be reconciled by adopting M, < O.53M0. 

No red edge was encountered in the low-log g instability 
domains. Due to our restriction of the computations on the 
long-period side, we could not follow overstable modes to tem- 
peratures below 125 000 K on the high-luminosity branch of 
the evolutionary tracks. At these longest periods that were com- 
puted, we did see no signs yet of a decline of the strength of the 
instabilities. 

The high-log g instability domains start at logTeff = 5.12 
for the O.53M0 model and at successively lower T^ffS for higher 
masses. If elemental diffusion was included, all instability do- 
mains showed a sharply defined red edge on the log Teff - pe- 
riod plane. The 0.59Mo model sequence was evolved once with 
diffusion at log T^ff < 5.0 and once without diffusion. The re- 
sulting instability boundaries are displayed as dot-dot-dashed 
lines in the left panel of Fig. |3 For log Tgff > 4.98 both in- 
stability domains are indistinguishable. The diffusive models 
encounter a red edge at logTeff ~ 4.97. The instability do- 
main of the diffusion-free sequence (thin dot-dot-dashed line), 
on the other hand, continues to grow towards lower T^ss with- 
out a sign of a weakening at log T^ff = 4.70, the location of the 
coolest model analyzed. 

The position of the red edge is not realistic; its location de- 
pends on the epoch when diffusion is switched on during the 
evolution computations. All that is important here is that dif- 
fusion causes a red edge; this point is further elaborated on in 
Sect. 10 

The maximal period- widths of the high-log g instability do- 
mains range from 290 to 480 s with no clear mass dependence. 
On the other hand, the magnitude of the shortest and the longest 
overstable periods are a function of mass. The 0.53Mq models 
have the longest-period overstable modes; the magnitude of the 
upper period boundary decrease with increasing mass. For the 
O.64M0 sequence, the longest overstable period is as short as 
500 s. On the other side, the shortest-period overstable modes 
have periods as low as 192 s for the O.64M0 sequence. This is 
shorter than the shortest observed periods in any GW Vir star. 

The most obvious problem is evidently the prototype of 
the class, PGl 159-035 itself. PGl 159-035 lies so close to the 
evolutionary knee of the model tracks that it falls mostly be- 
tween the low- and high-log g £ - \ instability domains. 
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Fig. 3. The boundaries of the instabihty domains on the log T^-sf - period domain. The left panel shows the dipole modes, the right 
panel contains the results for the quadrupole modes. The loci traced out by the four model sequences are plotted with different 
line types. For the 0.59Mo sequence, the results obtained on models with suppressed diffusion are shown in addition with a thin 
dot-dot-dashed line in the left panel. The boxes delineate the observational situation for dipole modes. Details are given in the 
text. 



The observational box on the log Teg - period plane overlaps 
only marginally with the instability domains of the O.SSM© se- 
quences. Such a low mass must be rejected, however, based on 
the computed period separations (see also Sect.l^. 

The right-hand panel of Fig. |3l shows the boundaries of 
the instability domains traced out by the quadrupole g modes. 
For the two most massive sequences, the instability topology 
looks the same as for the dipole modes. The main difference 
is that considerably shorter-period modes are overstable. For 
the 0.55 and 0.53Mq sequences, the instability domains are un- 
interrupted, they extend around the knee from the low- to the 
high-log g domain. The effective temperatures of the red edges 
of the dipole and quadrupole instability domains coincide. The 
growth rates of the overstable t - 2 modes reach comparable 
magnitudes as those of the H = I modes. 

For comparison with spectroscopically calibrated PG1159 
stars, Fig.|3plots the evolutionary tracks and the overstable re- 
gions on the log T^ff - logg plane. Observed GW Vir stars and 
the associated error boxes are superimposed. PNNVs are indi- 
cated by encircled dots. The loci of the blue edges of ^ = 1 and 
2 modes are given as dotted and dash-dotted lines, respectively. 
The evolutionary phases harboring overstable dipole g modes 
are marked with a heavier line on top of the evolutionary track. 
The general agreement between observations and models is 
reasonably good. Note again, the red edge position is not trust- 
worthy. Important is, however, the existence of a red edge. 

Finally, Fig. |5] shows the evolutionary tracks of the model 
stars on the HR plane; thereon, the blue edges can be approxi- 
mated by straight lines. The dipole-mode blue edge obeys the 
relation 

logL/Lo = 9.289 ■ log Teff - 45.186 . (1) 

The dash-dotted line in Fig. |5] delineating the quadrupole- 
mode blue edge, is less steep than the ^ = 1 blue edge and 



follows the relation 

logL/Lo = 9.0372 • log T^s - 44.158 . (2) 

The above parameterizations of the blue edges imply a mass 
dependence of the GW Vir instability domains as it is already 
evident in Fig. |3l In Sect. 0] we return to the observation that 
the orientation of the blue edge on the HR plane differs from 
that of the classical instability strip, in particular the slope as 
the opposite sign. 

3.1. The excitation meciianism 

The g modes that are found overstable in our models are all 
driven by the /c-mechanism induced by the opacity bump peak- 
ing around log T - 6.2. Partial ionization of K-shell electrons 
in C and O cause the peak. As the matter of the stellar en- 
velopes also contains heavy elements before the onset of diffu- 
sion, the high-temperature Z-bump enhances the opacity peak 
around log T - 6. Hence, the pulsation driving is the same as 
in all hitherto stability studies of PGl 159 variables. 

The abundance profile shown in Fig.|2lis representative of 
all models before the occurrence of diffusion-induced helium 
gradients in the outer envelopes (see Fig.|6jl. For the diffusion- 
less models, the abundances remain on the level as at the left 
border in Fig.|2]out to the photosphere. Note that in Fig.|2] the 
temperature reached already 1.6 ■ 10^ K at log(l -m/Mt) = -5. 
Hence, all model stars have at least 30.6% of helium (in mass 
units) in the driving region. In particular, all oscillating PGl 159 
models with log Teff > 5.0 (on the white-dwarf cooling tracks) 
have no composition gradient between photosphere and the 
pulsations' driving region. Ther efore, the ex c itation physics re - 
mains the same as discussed by Saio (' 1 996*). iGautschvl l l 1 997l) . 
and more recently bv .Ouirion et al..(.2004) . 
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Fig. 4. The GW Vir instability domains on the log TgSf - log g 
plane. The dots represent observed variable PG1159 stars. 
Encircled dots indicate the pulsators that are embedded in 
a planetary nebula. The boxes around the observational data 
mark the canonical uncertainties. The dotted line indicates 
the dipole-mode and the dash-dotted line the quadrupole blue 
edge. Along the evolutionary tracks, the epochs with overstable 
g modes are plotted with heavy lines. In the lower right of the 
figure, the two grey dashed lines mark the loci of constant pe- 
riod separation (21s and 22 s) as derived from the nonadiabatic 
stability computations. 



Fig. 5. The loci of all model sequences in the PGl 159 re- 
gion of the HR plane. The heavy lines delineate the do- 
mains of unstable f = I g modes. The blue edges for dipole 
and quadrupole modes, respectively, are indicated with broken 
lines. For logTeff > 5.1 at high luminosities of the 0.59Mo 
models, no instability is indicated; most probably, however, 
overstable ^-modes exist there. However, the computations 
were stopped because the shortest unstable periods exceeded 
1800 s. 



3.2. The red edge 

As mentioned before, all model sequences were evolved with 
switching on elemental diff'usion after they passed log T^ff = 
5.0 on the white-dwarf cooling track. The choice of the epoch 
of the switch-on of diff'usion was arbitrary. If diff'usion alone, 
without competing counteracting physical transport processes, 
is considered, sedimentation is very eff'ective. The top panel 
of Fig. |6l shows the helium profiles close to the g-mode driv- 
ing region at several epochs of the O.55M0 model sequence. 
It took the star to evolve only from logTgff = 5.0 to 4.951 
for the superficial He abundance to rise from 0.306 to 0.77. 
The superficial flattening of the He profile is an artifact of the 
modeling procedure: the diff'usion equations are solved only in 
the Henyey part of the stellar interior but not in the stellar en- 
velopes that are computed in full equilibrium. Hence our results 
probably underestimate the diff'usion effects. 

As mentioned in Sect. |3 we computed also the stability 
properties of 0.59Mq model stars devoid of diffusion. The cor- 
responding stability analyses were hence performed on models 
with a constant He mass abundance of 0.306 from the surface to 
the nuclear-active deep interior Neglecting difiiision altogther 
for the whole evolutionary history, no damping of the g-mode 
pulsations was found along the white-dwarf cooling track. 



On the other hand, if diffusion was accounted for, the 
g mode instabiUties literally starved when the He abundance 
exceeded a critical level. The microphysics of elemental trans- 
port is such that the He that is going to float on top of the 
heavier elements starts to separate from these other elements 
at a temperature depth of about 10^ K. The steepest abundance 
gradients are found between log T - 6.7 and 6.3, this region 
also dominates the damping and driving of the g-modes for the 
PGl 159 variables. Hence, it is understandable that superseding 
the main driving elements of that region with "pulsationally 
inert" stellar material reduces the efficiency of the pulsational 
driving. Looking at the differential work curves in the lower 
panel of Fig. |6l shows indeed that the curves are essentially 
self-similar. The higher the helium abundance in the driving 
region the less efficient is the driving contribution to the pulsa- 
tions. Finally, with Y= 0.77 in the superficial layers in Fig.|5J 
the driving cannot compensate the damping anymore and the 
g modes turn stable. The critical level of helium above the driv- 
ing region for the red edge varies between 0.77 and about 0.82 
for the different masses of the model stars considered here. 

The differential-work curves are very similar in form as the 
He abundance rises towards the surface; the compositional gra- 
dient does not distort the eigenfunctions (the spatial kinetic 
energy distribution is essentially unaffected by the diffusion 
process). Therefore, one can say that diffusion lets the pulsa- 
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Fig. 6. Top panel: Spatial variation of the He profile (Y) as a 
function of depth, measured in log T, for four 0.55Mo model 
stars as the red edge of the g-mode instability is approached. 
Bottom panel: Differential work, w = dW/dr, of the £ = l,k = 
19 mode for the same models. The eigenmode belonging to the 
thinnest hne is pulsationally stable. 



tions "starve". In the present computations, diffusion destroys 
the pulsations rather quickly after sedimentation is "switched 
on" at Jeff = 100 000 K. The red edges of the different mass se- 
quences were then encountered between 70 800 and 89 100 K. 
Approaching the red edge, the long-period modes died away 
earlier than the short-period ones. Therefore, close to the red 
edge, the mean period of the range of excited modes 



4. Discussion 

The splitting of the PGl 159 pulsation domains (see Fig. [Sjl 
into two regions - a long-period, low-log g domain and a short- 
pe riod, high-log pa rt - agrees with the results of lSaid ill 9961) 
and'Gautschv" (1997^. In the current computations, no continu- 
ous instability domain connecting low- and high-log g regions 
was found for dipole modes. In the two older studies just men- 
tioned, the 0.56 and 0.57Mo sequences had continuous instabil- 
ity domains, however The reason is connected with the differ- 
ent loci of the evolutionary tracks on the Teff - logg plane, the 
present models reach higher log g values earlier in their evolu- 
tion than in the old simplified ones of Saio and Gautschy. 



Regarding the comparison with observations, there is es- 
sentially one problem case: PGl 159 - 0354 itself. As men- 
tioned before, on the Teff - period plane, the observational box 
overlaps only slightly with the computed instability domains of 
the O.53M0 model sequence. Comparing the associated com- 
puted mean period-separation (Allf^i - 22.46 s) with the ob- 
served one (21.62 s) makes it clear that the O.53M0 model is 
no acceptable solution. With regard to the period-separation, 
the best agreement was found with the 0.55Mo model having 
Anf=i = 21.50 s at Teff = 138 000 K. This, however, means 
that the t - \ blue edge is at least 6000 K to cool there to 
be compatible with the observed situation. Our preferred mass 
fo r PGl 159 - 0354 i s clos er to the spectroscopic calibration 
of iDreizler & HebeJ (Il998l) (0.54Mo) than to the O.59M0 of 
iKawaler & BradlevI \\99m . The latter value is frequently re- 
garded as kind of a canonical mass for PGl 159 - 0354; it 
should, though, be kept in mind that the asteroseismic calibra- 
tions rely on fits to some evolutionary models with all their 
shortcomings and uncertainties. 

According to the current models, regions of exclusively 
t - 2 overstable modes were encountered for lower-mass (i.e. 
about 0.53 < M, < O.57M0) PWDs when they pass along the 
low-luminosity side of the evolutionary knee. 

As mentioned hitherto for DBVs jGautschv & Althau3 

2002), also for the pulsating PGl 159 models we found the 
shortest-period overstable modes to be quadrupole ones and 
the longest-period overstable mod es to be dipole ones . Even 
if the latest RXJ21 17 observations ( Vauclair et alJ20()2 ^ do not 
clearly reveal quadrupole modes, the most likely € - 2 can- 
didates in the data are at the short-period bound ary of the ob- 
served set. Also for PGl 159 ( Wing et et alJl99ll) . the positively 
identified ( - 2 mode has a short period at 425 s. However, one 
multiplet, the one associated with longest-period mode, is pos- 
sibly due to an /" -2 mode; if true, this would clearly contradict 
our computations. Another potential problem case is PG0122; 
its 465 s mode - th e third-longest period mode reported in 
O'Brien et al. ( 1996) - is suggested to be a quadrupole mode. 
The identification was not without debate, however. 

With regard to mode excitation, we saw that neither the use 
of overly simplified full stellar models nor the use of envelope 
models only are the reason for overstable g modes in envelopes 
with helium mass abundances of about 0.3. Additionally, we 
have now three different numerical methods (finite-difference 
relaxation, finite-element solver, and a nonlinear shooting 
method) that a gree on the ;?-mode excitation. The disagree- 
ing results (e.g.lstarrfield et al.ll983tlBradlev & Dziembowskil 
Il996t Ic^HooT were obtained with finite-difference relax- 
ation methods which resemble the method of Saio ( 1996). We 
hypothesize therefore that the origin of the disagreement is 
linked to particularities of the computational realizations. We 
found the previously proclaimed "helium poisoning" to be- 
come effective only at higher mass abundances such as they 
are encountered when elemental diffusion was at work. If the 
relative mass abundance of helium exceeded about 0.77, the 
g modes were found to "starve" so that a red edge occurred. 

With respect to pulsational driving, the only qualitative dif- 
ference to the results of .S aia (,1996.) and .Gautsch v. 11997.) is the 
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absence in the current model stars of e-driven instabilities with 
periods of around 100 s. 

How can the bimodal period distribution on the T^s - pe- 
riod plane (cf. Fig.|3} be understood? At a chosen Teg, the peri- 
ods differ by at least a factor of three, whereas the luminosities 
differ by factor of 15 or more. In contrast to the expectations, 
the short periods are excited at low and long periods at high lu- 
minosities. To understand this, we resort to the classical time- 
scale argument (e.g. in Cox & Giuli 1968) 



/ 

JAi 



CvTAmlU = 0{ni2n). 



(3) 



The integral extends over the mass of the envelope, AM, over- 
lying driving region of the pulsations. Relation Q defines the 
pulsation period, 11, for which the thermo-mechanical coupling 
in the stellar material is optimal for pulsational driving. We 
compared two 0.59Mq models at logTeff - 5.06, one at high 
and the other at low luminosity. For both models, we evalu- 
ate the left-hand side in eq. Q. We approximate the integral, 
as usual, by adopting the values at the driving region for the 
quantities cy and T (indicated by a tilde) and write it hence 
as cv T AM. Based on the full stability computations, we as- 
sumed the driving to occur at around log t = 6.3 in all cases. 
Therefore, we compute the ratios of cV,AM, and L,; all three 
quantities are considerably larger than unity as the evolution 
proceeds past the knee (always measured at the same Teff). It 
turns out that the product in the numerator of Eq. Q over- 
compensates the change in L» in the denominator. Therefore, 
the high-luminosity phase of evolution prefers longer periods 
for pulsation than the low-luminosity one. The overcompensa- 
tion is mainly due to the considerable decrease of the envelope 
mass, overlying the driving region, as the stars evolve towards 
to white-dwarf cooling track. At log T^ff - 5.06 for example, 
AM at high L* is 630 times larger than at low L,. The influence 
of this number is further enhanced by the fact that radiation 
pressure constitutes a considerable fraction of the total pres- 
sure, P, in the high-luminosity stage. Since cy oc Ijff-, with 
P — \ - PailP, cy can rise considerably when radiation pres- 
sure gains importance. Again, for the models at log Teff - 5.06 
we found yS =s 0.4 at the high luminosity and ^ a; 0.9 at the 
low-luminosity epoch. 

To understand the slope/orientation of the blue edge in the 
post-knee domain in the HR diagram, we follow the reason- 
ing of Cox & Hansen ( 1979). We approximate eq. the same 
way as to understand the period bimodality. We approximate 
the hydrostatic equilibrium as 



\tiPR'^ 
GM 



AM. 



The Stefan-Boltzmann relation allows to rewrite eq. Q as 



cvTPL 
n T\. M 

eri 



: const. 



The envelope P - T structure of PWDs is decently approxi- 
mated by a radiative-zero one using a Kramers-type opacity. 



Eliminate P in the last equation with the chosen P ■ 
to get 

nTl^ ■ VM 



const. 



T relation 



(4) 



For the period H, we adopt the relation D oc T,„-R from 
Winaet et al. '('1983'). Using a very simple white-dwarf model 
from Kippenhahn & Weiaert ( 1994) (their section 35.3): oc 
(L/M)^^^ to estimate the maximum temperature, T„„ leads to 

UccR ( — ] ' . 

\mI 

Eliminating R again via Stefan-Boltzmann's law gives 



CyT'-^ LI 

Mt* 



'7^6 

■'efF 



const. 



(5) 



Assume the temperature at the top of the driving region, f, to 
remain constant within the instability domain of the PGl 159 
variables and neglect the mass dependence of the derived re- 
lation, as the range of stellar masses assigned to PGl 159 vari- 
ables is small. From fits to the evolutionary computations we 
derive dlncV/dlnLlj-^^ ^ 0.5 and we neglect any additional 
Teff dependence. All things considered, we finally arrive at: 



dlnL 



28 

y^9.3. 



(6) 

dlnTeff 3 

According to eq. Q, the slope of the blue edge as obtained 
from the fit to the full nonadiabatic computations is 9.3 for 
{ = I and 9.0 for £ - 2. Regarding the simplicity of the model 
and the coarseness of the approximations, the result is aston- 
ishing. We attribute the dominant contribution to the slope of 
the blue edge to the term d In cy/d In L. Hence, it is the growing 
importance of radiation pressure as the luminosity increases in 
the PGl 159 envelopes that tilts the blue edge to the left on the 
HR- and the log T^s - \ogg diagrams. 

We found elemental diffusion to quench g-mode instability 
on the high-log^ branch of evolution when the helium abun- 
dance exceeds a critical level in the driving region. A helium 
mass abundance of about 0.77 was found to be sufficient for 
the overstable s r nodes to starve. Hence, as conjectured by 
IPreizler & H eber (1998), a depletion of C and O, which can 
be regarded complementarily as an enhancement of He, due to 
sedimentation in the driving region enforces a red edge. 

Diffusion is a very efficient process to turn PGl 159 stars 
into DO stars. If only diffusion - without counteracting pro- 
cesses - is included in stellar modeling, DO stars would show 
up at too high effective temperatures to be compatible with 
observations. In the current modeling, diffusion was switched 
on ad hoc at 100000 K on the white-dwarf cooling track. 
Therefore, the location of the red edge of the PGl 159 variables 
in ffie various diagrams shown in this paper is not trustworthy. 

The quantitative non-equilibrium computations of 
lUnalaub & Bues' (200&) to model the counteracting processes 
of diffusion and mass loss suggest that a unique red edge for 
the PGl 159 variables as a class cannot be expected. As com- 
positions and rotational velocities (which tend to homogenize 
composition gradients) vary from star to star, also the epoch 
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(and hence the corresponding Tf.fi, L, or g) will vary when 
the critical contamination of the driving region is reached. 
Therefore, it is not inconceivable that close to the red edge of 
the PG1159 variables, these are intermixed with stable DO 
stars and even with non-pulsating PGl 159 stars. The closer 
one gets to the red edge, the shorter should the mean period 
of the excited modes be; this occurs since, according to our 
computations, the long-period modes starve earlier than the 
short-period ones. 

The position of the theoretical instability domains is only 
part of the game to understand the GW Vir pulsators. Pertinent 
information is obtained from period-spacing data. We used 
the results from ( - \ computations to compare them with 
the observed period separations for three stars. In the case of 
RXJ21 17 H- 3412, only tracks of 0.64 and 0.59Mq models pass 
on the log Teff - log g plane within the observational error box. 
Inspecting the period separations revealed that the 0.64 mod- 
els have too small a mean period separation. Even the O.59M0 
objects with the arithmetic-mean period separation Anf=i « 
20.0 s were below t he mean observed separation of 21.62 s 
llVauclair et al.l20()2l) . Hence, only a mass below O.59M0 can fit 
the observation. However, the corresponding current evolution- 
ary tracks pass outside the error box. PGl 159 - 035' s observed 
AII^=i = 21.62 s period spacing JWinget etal.ll 99 ih is best re- 
produced with our O.55M0 models. As mentioned before, the 
problem with the O.55M0 model stars is that none of the in- 
volved dipole modes is pulsationally overstable in the vicinity 
of the position of PGl 159 - 035 on the log Teff - \ogg plane. 
Again, evolutionary models that arrive at the lower side of the 
knee with higher luminosity than what we get currently might 
resolve the dilemma. Finally in PG2 1 3 1 + 066, aH^ = 2 1 .65 s 
is observed (.Kawaler et al. 1995 ). From our current modeling, 
we conclude that the best fitting stellar mass lies between 0.55 
and O.59M0; the observed pulsation modes fall then well within 
the computed range of overstable dipole g modes (cf. Fig.|3}- 

Figure0]has inscribed lines of constant period separations. 
The loci of 21 s and 22 s, interpolated from the nonadiabatic 
computations, are shown. Note that the observed period sep- 
aratio ns for RXJ2117 , PGl 159 - 035, PG2131, and (accord- 
ing to IO'Rrienll200nl) PG0122 all he between 21 and 22 s. If 
PG1707 + 427's period separation of 18.37 s jFontaine et alJ 
I199I,) holds then the apparent confinement to 2 1 - 22 s of the 
other four GW Vir variables is an observational bias. However, 
in case the period separations of GW Vir stars are indeed nar- 
rowly confined, then they populate only a fraction of the insta- 
bility region as computed by linear theory. It appears worth- 
while to point out that the loci of constant period separation 
depend also on the trapping amplitude of the modes. Our mean 
period separations were computed as simple arithmetic means 
over many successive radial orders. The larger the trapping 
amplitudes are, i.e. the sharper the composition gradients in 
the stellar interior, the bigger the difference between AH and 
the asymptotic value. If mass-dependent diffusion efficiencies 
are capable to increase the distance between the, say 21 and 
22 s, period- separation lines in Fig. |3 remains to be seen. 
On the observational side, the number of asteroseismically 
studied PGl 159 is so low, that every additional well-resolved 



oscillation-frequency spectrum that becomes available makes 
an important contribution to solve this conundrum. 

In addition to the period-separation information, we can 
also compute the lengths of the trapping cycles. The two com- 
parisons with observations that are currently possible revealed 
the modeled trapping -cycle lengths to be always shorter than 
the observed ones: RXJ2117 -1- 3412, 83.9 s observed and 
56 - 60 s computed. For PGl 159 - 035: 80.5 s observed and 
52.6 - 68. 4 s computed. Based o n mode-trapping theory as de- 
veloped in lBrassard et alJ lll992h . the current model stars have 
seemingly too thick helium layers to be compatible with obser- 
vations. As there is ample evidence for mass loss before and 
during the PGl 159 episodes of evolution ( Werner. .200 L) . we 
are positive that the mass of the helium-rich model envelopes 
can be sufficiently reduced via a hot stellar wind to improve the 
agreement with observations. 

5. Conclusions 

The exc itation phys ics a s scraped out fro m simplified model 
stars by Saio (1996) and"Gautschvl lll997h is corroborated us- 
ing now complex stellar-evolution models that went through 
a very late thermal flash. In particular, we found again that 
modes can be excited with considerable He abundance in the 
driving region, supporting the view that abundance gradients 
between the spectroscopically accessible stellar surface and the 
driving region of g modes do not have to be invoke d to explain 
PGl 159-type oscillations (e.g. Bradlev & Dziembowski 1993) . 

A blue edge was encountered for dipole and quadrupole 
modes; its slope is a manifestation of a mass dependence 
of the pulsational instability of PGl 159 stars. On the high- 
luminosity branch of the evolutionary track, the hottest vari- 
ables are the most massive ones. The situation reverses on the 
low-luminosity branch (below the evolutionary knee), where 
the bluest variables are the least massive ones. 

Model-star sequences including diffusion led to a red edge 
of the PGl 15 9 instability domain . Much in agreement with the 
predictions of iDreizler & HebeJ (0.998), the oscillations of the 
PGl 159 stars starve once the He abundance in the driving re- 
gion exceeds a critical magnitude. It is therefore likely that the 
red edge is close to the transition from PGl 159 to hot DO- 
or possibly DA-type white dwarfs. It is furthermore conceiv- 
able that the red edge of the GW Vir stars is not sharply de- 
fined but rather frayed. The counteracting effects of rotation, 
stellar wind, radiative levitation and gravitational settling de- 
termine for each star individually when the g-mode pulsations 
die out. A significant step ahead would be parameter studies of 
the competition of mass-loss and/or diffusion on the efficiency 
of pulsational driving. Currently, diffusion is still switched on 
ad-hoc. 

In contrast to the models of 'Saio' ('1996') and 'Oautschvl 
(Il997t) . the stability analyses of the current stellar model se- 
quences did not reveal any overstable e-mechanism driven 
g modes at short periods. 

We gathered some evidence, when comparing the theoret- 
ical with the observed period- separation data, that the current 
evolutionary tracks evolve at too low luminosity on the HR dia- 
gram and hence on the log Teff-log g plane. Furthermore, based 
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on the magnitudes of the trapping-cycle lengths, the helium- 
rich envelopes of the current PG1159 models seem to be too 
massive to agree with observations. It appears therefore worth- 
while to look again into the microphysics and the numerics of 
evolutionary computations to try to resolve the persisting dis- 
crepancies. 
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